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The inner boundary of a black hole accretion disk is often set to the marginally stable circular 
orbit (or the innermost stable circular orbit, ISCO) around the black hole. It is important for the 
theories of black hole accretion disks and their applications to astrophysical black hole systems. 
Traditionally, the marginally stable circular orbit is obtained by considering the equatorial motion 
of a test particle around a black hole. However, in reality the accretion flow around black holes 
consists of fluid, in which the pressure often plays an important role. Here we consider the influence 
of fluid pressure on the location of marginally stable circular orbit around black holes. It is found 
that when the temperature of the fluid is so low that the thermal energy of a particle is much 
smaller than its rest energy, the location of marginally stable circular orbit is almost the same as 
that in the test particle case. However, we demonstrate that in some special cases the marginally 
stable circular orbit can be different when the fluid pressure is large and the thermal energy becomes 
non-negligible comparing with the rest energy. We present our results for both the cases of non- 
spinning and spinning black holes. The influences of our results on the black hole spin parameter 
measurement in X-ray binaries and the energy release efficiency of accretion flows around black holes 
are discussed. 
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I. INTRODUCTION 


In the spacetime around a black hole, the circular mo- 
tion of a test particle is not always possible according to 
general relativity. There is a smallest radius on which 
the circular motion of a test particle is marginally sta- 
ble p. This radius is called marginally stable circular 
orbit. With a small perturbation, the circularly moving 
particle at this radius will then plunge into the black hole 
freely. The properties of the transition from inspiral to 
plunge depend on the mass ratio 7 of the particle and 
the black hole [2 The transition would be less_grad- 
ual (or more ”abrupt”) with smaller 7 (See also Bid 
This property is more relevant in the context of the in- 
spiral in a binary black hole system, when the two black 
holes are just about to merge. When the self-force of 
the inspiralling particle is considered, the location of the 
marginally stable orbit is also modified [5]. 

The marginally stable circular orbit is also important 
in the estimate of energy release of the accretion to a 
black hole. By calculating the binding energy of the cir- 
cular motion on this radius |6], one can estimate how 
much of the total energy can be released during the ac- 
cretion process. For the accretion flow around a black 
hole, the energy release efficiency rangs from 5.6% (for 
non-spinning or Schwarzschild black holes) to 42% (for 
extreme Kerr black holes). The location of the marginally 
stable circular orbit and the energy release efficiency are 
closely related, both depending on the spin parameter of 
the black hole [7]. 

In some applications of accretion disk model in astro- 
physical observations, the location of marginally stable 
circular orbit is also crucial. Since the marginally stable 
circular orbit depends on the spin of the central black 
hole, it is possible to measure the spin of a black hole 


if the corresponding marginally stable orbit can be mea- 
sured. By fitting the observed soft state spectra of black 
hole X-ray binaries, which are assumed to come from a 
thin disk with its inner boundary at the marginally sta- 
ble circular orbit, one can derive the location of this orbit 
and estimate the spin parameters of the black holes in X- 
ray binaries |8]. 

In the theory of accretion flow, the marginally stable 
circular orbit is also important. Due to its transitional 
nature from the inspiralling region to the plunging re- 
gion, it is often believed that an accretion disk is torque 
free at this radius. Although still being debated when 
considering magnetic fields Aid, the torque free condi- 
tion on the marginally stable circular orbit serves as an 
important inner boundary condition in many accretion 
disk models. 

The self-gravity of the particle is not relevant in the 
context of accretion disks, but the material in a real ac- 
cretion disk is fluid rather than test particles. When 
considering fluid, the particles within it are interacting 
with each other, and the pressure plays an important 
role in the dynamics. In the study of accretion disk the- 
ory, the marginally stable circular orbit for a test par- 
ticle is used as the inner boundary of an accretion disk 
in many models. More precisely, it has been shown by 
analytical theory of fluid tori around black holes that 
the inner boundary of an accretion disk lies between the 
marginally stable circular orbit and the marginally bound 
orbit, while both expressions are for test particles HHz. 
In this work, however, we focus on a different problem, 
that is, we try to investigate the marginally stable cir- 
cular orbit itself, with the influence from the pressure in 
the fluid. In section 2, we briefly mention the marginally 
stable circular orbit in the case of a test particle. Then 
we consider a thin disk consisting of perfect fluid around 


a black hole and discuss the marginally stable circular 
orbit in this case in section 3. In Section 4 we present 
a brief discussion on our results. In this paper, we set 
c= G= M = 1, with c, G, and M the speed of light, 
gravitational constant, and the mass of the black hole, 
respectively. We use (+ — — —) signature all through. 


II. MARGINALLY STABLE CIRCULAR ORBIT 
OF A TEST PARTICLE 


The metric of the spacetime outside a Kerr (spinning, 
non-charged) black hole is 
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For a steady and axis-symmetric spacetime (e.g. 
Schwarzschild, Kerr), there are two apparent constants 
of motion of a free particle with unit mass 


u=E, ug =—-L (4) 


Since the metric is block diagonal, the normalization con- 
dition of the 4-velocity can be written as 
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Consider the circular orbital motion on the equatorial 
plane (0 = 7/2, u? = d0 /dr = 0), the equation of motion 


can be written as 
2 
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In order to maintain a circular orbit, dr/dr and 
d?r/dr? should both vanish, which imply 


U, = 0, U, = 0, (8) 


where the prime denotes the derivative to r. The two 
constants of motion can be derived from the above two 
equations as 
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For prograde motion, a > 0; for retrograde motion, a < 0. 
Not all the circular orbits are stable. The stable ones 
should also fulfill U/’ < 0, and the marginal stable cir- 
cular orbit corresponds to U/’ = 0. Using the expression 
for E and L, we can derive the marginal stable orbit 


Tms = 3 + Z F [(3 - Z1)(3 + A +222), (11) 


where 
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The upper sign and the lower sign are for prograde and 
retrograde motion of the particle, respectively. Note that 
for a Schwarzschild black hole, a = 0, rms = 6. The ex- 
pression above has been widely adopted in various astro- 
physical studies of black hole systems. 


II. MARGINALLY STABLE CIRCULAR ORBIT 
OF A PERFECT FLUID DISK 


For perfect fluid, the energy-momentum tensor can be 
written as 


Ty = (p + E)u"uy + por, (12) 


where u”, p, and £ are the 4-velocity, pressure, and the 
energy density respectively. It can be proved that there 
are also two constants of motion (13): 
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and 


ug = -L, (14) 


where n is the number density. Note that these two con- 
stants of motion represent no longer the energy and angu- 
lar momentum. Use the normalization of the 4-velocity, 
it is easy to show that 
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where the effective potential 
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A. The Schwarzschild Black Hole Case 


We first consider the case of a thin fluid disk around 
a non-rotating black hole (with spin parameter a = 0). 
The effective potential becomes 


U, = (; + :) er — (r? aryl] — (rf = 2r°). (17) 


In this case, the energy density £ consists of two parts: 
One is the rest energy (proportional to the number den- 
sity n) and the other one is the thermal energy (assumed 
to be proportional to the pressure p). Therefore we have 
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E= Mn+ Ti” (18) 


where mo and y are the rest mass of the particle and the 
ratio of specific heat, both of which are constants. If the 
dependence of both number density and pressure on the 
radius are of power law forms, we can parameterize the 
function n/(p + €) as 


7 = 75, (19) 
pte 14+Cr-? 
where we have assumed p/p x T œ r~°, where p = mon 
is the density of the fluid. For the standard thin disk 
model [7,14], bis 3/8, 9/10, 3/4, in the inner, middle and 
outer regions, respectively. For the self-similar solution 
of an advection dominated accretion flow (ADAF, [15, 
[16]), b equals to 1. However, one should note that these 
dependence on radius r is only for the region where r >> 
Tms. When r ~ rms, b can be negative if the torque 
free condition at r ~ rms is applied, which may be more 
relevant in the case of real accretion flows around black 
holes. The constant C is always positive. 

Since the constant B can be absorbed into € and £, 
the effective potential can be rewritten as 
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In a normal fluid, the temperature is usually low and 
the thermal energy is much smaller than the rest energy, 
that is 


= pP MN. (21) 

y—1 
In another word, the constant C is very small (C < 1) 
in this case. Then there is not much influence from the 
pressure. The effective potential is approximately the 
same as that in the test particle case, which means the 
location of the marginally stable circular orbit is almost 
the same even in the fluid case. 

However, in some hot accretion flows around black 
holes, e.g. ADAF (15), [16], the thermal energy may be 
relevant. We can estimate the ratio of thermal energy to 
rest energy in this case as 
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where k, Ti, Mmi, Te, Me, c are the Boltzmann constant, 
ion temperature, ion mass, electron temperature, and 
speed of light, respectively. In this case, C is not very 
large (C < 1) but also non-negligible. We can expand 
the marginal stable circular orbit as 


Tms = Tms,0 + CArms.- (24) 


In order to calculate the correction term Arms, we have 
to know the first and second derivatives of the effective 
potential, Eq. (20). They are 
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Setting Eq. (20) and above two equations to 0 (U, = 
U! = U” = 0), we get an equation after eliminating E 
and L, 
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Keeping the terms of the order O(C) we can get 
Arms = brine o{(Tms,o — 3)[(7 — b)rms,o + (2b — 10)] 
—(Tms,0 = 2)(5rms,0 pe 12))} 
(28) 
For a Schwarzschild spacetime, r5,9 = 6, so we get 


Tms = Tms 0 + CArms =6+6 °Cb(24—12b). (29) 
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FIG. 1: The dependence of the marginally stable circular orbit 
Tms in Schwarzschild spacetime on the power-law index b in 
the case of C < 1. The solid line and dashed line are for 
C = 0.1 and C = 0.3, respectively. The dotted horizontal line 
represents rms = 6 in the test particle case. 


In the equation above, Arms(b) has its maximum 3.75 
at b = 0.413 and local minimum —0.18 at b = 2.7, respec- 
tively. So the maximum of the marginally stable circular 
orbit is 


T'ms,max = 6+ 3.75C, (30) 
corresponding to b = 0.413 and the local minimum is 


Tms,min = 6 — 0.18C, (31) 


corresponding to b = 2.7. When b > œ, rms — 6, and 
when b is negative, rm; < 6 and it decreases monotoni- 
cally with the decreasing of b (see Fig. I). Therefore, in 
this case when the thermal energy is smaller than the rest 
energy, the influence of the pressure on the marginally 
stable circular orbit is non-negligible, especially when 
b< 1.5. 
B. The Kerr Black Hole Case 


In Kerr spacetime, the effective potential is expressed 
by Eq. (16). When C is small, we can do the same 
expansion as in Schwarzschild spacetime. 

In Kerr spacetime, the effective potential, and its first 
and second order derivatives are 
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respectively. After some long but straightforward deduc- 
tions, we can get 


Armes = 


where ro and Eo are expressed as 


ro = 3 + Z2 F [(3 - Z1)(3 + Zı + 222)", 


and 
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respectively, where 
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Zi=1+4 (1-a) (1 +a} + (1 - a)®], 
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E} is the derivative of Eo with ro. E? is the correction 
term of E? with order O(C). 

The behavior of the correction term is similar to that of 
the Schwarzschild case. This can be seen in Fig. [2] As an 
example, Fig. B]shows the dependence of the marginally 
stable circular orbit on b and C in Kerr spacetime with 
a = 0.5. 
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FIG. 2: The dependence of the correction term Arms on the 
power-law index b and spin parameter a. Dotted line: a = 0; 
solid line: a = 0.1; dashed line: a = 0.5; dash-dotted line: 
a = 0.9. 


IV. DISCUSSION 


The marginally stable circular orbit can be treated as 
the inner boundary of an accretion disk around a black 
hole in some cases. This is important because it is cru- 
cial to study the structure of accretion disks. It is also 
important when calculating the radiation of a standard 
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FIG. 3: The dependence of the marginally stable circular orbit 
Tms in Kerr spacetime with a = 0.5 on the power-law index b 
in the case of C < 1. The solid line and dashed line are for 
C = 0.1 and C = 0.3, respectively. The dotted horizontal line 
represents Tms = 4.233 in the test particle case for a = 0.5. 


accretion disk, where it is usually assumed that there is 
no radiation from the region inside the marginally stable 
circular orbit. 


However, the widely accepted marginally stable circu- 
lar orbit is based on the test particle assumption. When 
the temperature in an accretion disk is not very high 
(namely the thermal energy is much smaller than rest 
energy), the expression of the marginally stable circular 
orbit for a test particle is accurate enough. But when the 
temperature is so high that the thermal energy is non- 
negligible comparing to the rest energy, pressure may in- 
troduce some corrections to the location of the marginally 
stable circular orbit. 


If the torque free boundary condition is applied at the 
inner boundary of accretion disks, the temperature near 
the marginally stable circular orbit would decrease with 
the decreasing of radius, which is different from the region 
far from the inner boundary. Assuming the dependence 
of temperature on radius near the marginally stable cir- 
cular orbit is a power-law T œ r~°, then b < 0. As can 
be seen from Fig. [2] the correction term changes signif- 
icantly with b when b is negative. As an example, for 
the Schwarzschild case, when C = 0.2 and b = —0.2, the 
marginally stable circular orbit is rms = 5. 


The correction to the marginally stable circular or- 
bit not only influence the inner boundary of accretion 
disks, but also affect the energy release of the accretion 
flow around black holes. In Schwarzschild spacetime, the 


binding energy E (see Eq. D) at a certain radius r is 


r—2 

E= Rp ay (38) 
which equals to 0.9428 at r = rms = 6, so the correspond- 
ing efficiency is 1 — 0.9428 = 5.72%. If the marginally 
stable circular orbit is different, the efficiency can be dif- 
ferent, e.g., at r = 5 the efficiency is 5.13%. In Kerr 
spacetime, the energy E of a test particle co-rotating in 
a circular orbit at r is described by Eq. (9) The depen- 
dence of the efficiency 1— E on radius r is shown in Fig. 
Note that it is not monotonic, when the radius r < rms, 
the efficiency is also lower than that at the marginally 
stable circular orbit. 
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FIG. 4: The dependence of the efficiency 1 — E on radius r. 
The solid line, dashed line anddash-dotted line correspond to 
the cases of a =0, 0.5 and 0.9, respectively. 


However, not only the true location of inner boundary 
of accretion disks is different from marginally stable cir- 
cular orbit as shown by the analytic models ime) but 
as also mentioned, in a disk where the advection is 
important, it is not proper to use the marginally stable 
circular orbit as its inner boundary to calculate the ra- 
diation of the disk. The reason is that in an advection 
dominated flow, there is no equilibrium” of the circular 
motion, and the condition based on which the marginally 
stable circular orbit is derived is not fulfilled. But the 
widely used technique to measure black hole spin by fit- 
ting thermal soft state spectra ig, is all right because 
in a standard thin disk, the marginally stable orbit is al- 
most the same as for the test particle case, which is used 
for the measurement of the black hole spin parameters in 
X-ray binaries. 

One thing should be noted is that our treatment is 
based on a perfect fluid disk. However, a real accretion 
disk consists of viscous fluid. Another point that may 
worth further consideration is the influence of pressure 
on the ”abruptness” of the transition from inspiral to 
plunge regions. Because this property will heavily affect 
the inner boundary condition of accretion disks, if the 
transition is somehow less ” abrupt”, the fluid disk may 
continue further in beyond the marginally stable circular 


orbit. However, this is out of the scope of this work. We 
hope to tackle this problem in a future work. 
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